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We discuss the scaling of characteristic lengths in diffusion limited aggregation (DLA) clusters in 
light of recent developments using conformal maps. We are led to the conjecture that the apparently 
anomalous scaling of lengths is due to one slow crossover. This is supported by an analytical 
argument for the scaling of the penetration depth of newly arrived random walkers, and by numerical 
evidence on the Laurent coefficients which uniquely determine each cluster. We find a single crossover 
exponent of —0.3 for all the characteristic lengths in DLA. This gives a hint about the structure of 
the renormalization group for this problem. 

PACS numbers; 64.60.Ak, 61.43. Hv 



Diffusion-limited aggregation (DLA) is a model for 
growth that was introduced in 1981 by Witten and 
Sander [|l| and has been the subject of intensive scrutiny 
ever since. Remarkable progress has been made in 
numerical treatments of the model and in applications 
to various physical problems. Nevertheless, a fundamen- 
tal understanding of the most striking features of DLA 
clusters is still lacking. A very promising recent devel- 
opment is the method of iterated conformal maps intro- 
duced by Hastings and Levitov |3| (HL) and pursued by 
Davidovitch et al. Q. Using the method of conformal 
maps we will try to clear up one of the major difficul- 
ties in the field, namely the large corrections to scaling 
which lead to slow crossovers. This discussion has a very 
significant outcome: there has been serious doubt about 
whether DLA clusters are really fractal, based on the 
existence of characteristic lengths apparently scaling dif- 
ferently from the overall cluster radius. We show that 
these lengths cross over to scale in the same way as the 
radius, and suggest that all of the various deviations from 
exact fractal scaling have the same source, with the same 
correction to scaling exponent. There is, in our view, no 
longer any reason to doubt that DLA clusters are fractal. 

In the DLA model a nucleation center is fixed at the 
origin of coordinates and random walkers are released 
from outside. When a walker comes in contact with the 
cluster it sticks. The process continues until N walkers 
have attached; in modern work A^'s of 10^ are reasonably 
easy to attain. Provided lattice anisotropy is avoided, the 
clusters so produced seem to be self-similar with fractal 
dimension D w 1.71 in two dimensions and the proba- 
bility of growth at a point on the surface appears to be 
multifractal. 

The conformal map method uses the Laplacian growth 
version of DLA. That is, we take the cluster to be 
a grounded conductor in the z plane, with a probabil- 
ity to grow at a point on its surface proportional to the 
charge there: | Vl^|, where V is the potential with bound- 
ary conditions of unit flux at infinity and = on the 



surface. If we construct a complex potential such that 
Re[^{z)] = V, and define u; = e* = Z-^{z), then Z{w) 
is a conformal map which takes the exterior of the unit 
circle in the w plane to the exterior of the DLA cluster in 
the z plane, and is linear, Z ^ roW, for large \w\. Thus 
V \n{\z\/ro) at large \z\. The Laplace radius, Tq, is 
the radius of the grounded disk with the same capaci- 
tance (with respect to any distant reference point) as the 
cluster The growth probability is 1/\Z'\. Intervals dO 
on the unit circle in the w plane correspond to intervals 
of arc length ds with equal growth probability in the z 
plane. 

To characterize the map we write: 

Z{w)^roW + Y,Ak/w'' (1) 

k=l 

The Laurent coefficients, Ak, are a useful parameteriza- 
tion of the map. In what follows we will assume that the 
Laurent expansion of Z has no constant term. This cor- 
responds to moving the DLA cluster so that the center 
of charge is at the origin. 

In the HL method the Ak are produced directly. 
However, this method employs certain approximations, 
and also numerically difficult for large N because it is of 
order iV^, and as a practical matter is limited to N ^ 10''. 
We observe that there is a way around these problems. 
Suppose we produce a DLA cluster by the conventional 
method with random walkers (which is much faster) and 
then freeze it at some N . Then by recording where M 
random walkers would attach to the cluster we have a set 
of points 2„i- These are at angle 9m ~ 2-K{m/M) in the w 
plane, since we are sampling the charge, and equal incre- 
ments of charge correspond to equal increments of 9. The 
Ak are the Fourier coefficients of the function z{9,n)- We 
should note that everything in this paper results from the 
existence of the conformal map, and the behavior of the 
Ak, not on the detailed method of generation proposed 
by HL. 

Numerically we have found, using either the HL 
method 01 or the conventional method described above 
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(cf. Fig. |), that To cx N'^'^ and that for all but the first 
few fc's, (|Afc|2) cx N'^/^. For fc < 4 the At appear to 
scale with a smaller power of A^, a behavior that we will 
interpret below as a crossover. 

We will be concerned first with penetration depths of 
the random walkers, i.e. how far from the origin they 
land. This has been a matter of considerable interest for 
some time because the width of the distribution of this 
quantity ("the width of the growth zone") seems to in- 
crease more slowly with N than the mean radius . This 
is a disturbing observation since a real fractal should have 
no length scale other than its overall size. Some authors 
1^ have maintained that what was observed was a slow 
crossover and that asymptotically all lengths scale to- 
gether, but others [H have given evidence for an 'infinite 
drift scenario' in which DLA is not a fractal at all in the 
asymptotic limit. The structure of conformal map theory 
gives us a very elegant way to discuss these matters. 

We can define a penetration depth as follows: suppose 
we follow a field line from large |z| to the surface of the 
cluster. The displacement of the endpoint of the line 
from where it would terminate on the equivalent disk 
may be written as (r|| -I- ir±)w = Z{w) — tqw, where 
r|| is the radial displacement and r±_ is the transverse 
displacement, see Figure |[ The spread of r|| defines a 
penetration depth ^||: 



We also consider the analogous quantity ^_l: 



Using Eq. (|l|) we have 



(2) 



(3) 



(4) 



This exact result has a useful geometric interpretation: 
penetration of a DLA cluster corresponds to landing in- 
side one of the "fjords". Clearly r± is the azimuthal 
deviation of a walker from a straight path, and r|| the 
penetration. Since they are the same, on average, the 
only way for DLA to have an anomalous scaling for the 
penetration depth is for the angular width of the fjords 
to be smaller and smaller as the cluster grows. 

Now we can discuss the scaling of ^|| , the quantity of 
most interest, by analyzing the Ak and using Eq. (Q). 
We have remarked above that the numerical evidence is 
that the most of the \Ak\ scale as the radius. This is, in 
some sense obvious, since there is an elementary formula 
that relates the area, Sn, of the image of the circle to 
the Laurent coefficients: 



S 



N 



TrJ2k\Ak 



(5) 



The right hand side scales as N'^/^ , but the area only as 
N . The leading behavior of the sum must also be N^^-^ . 
The obvious way for this occur is if each coefficient has 
this scaling, which would lead to cx r^. We will now 
show that this expectation is correct. 

If we start with all the Ak small, then in the course 
of growth, they will not remain small. This is because 
Laplacian growth is subject to the MuUins-Sekerka in- 
stability ijT^. In fact the classic calculation of Ref. px| ] 
was concerned with exactly this behavior: the Fourier 
coefficients of the wrinkling of the surface grow expo- 
nentially, with the growth rate proportional to k. The 
slowest growth is associated with the smallest k. How- 
ever, we are concerned here with the limiting behavior of 
Ak, far out of the linear regime. 

We can get some insight by making the following esti- 
mate for large k: 

J d6'id6l2Z(e*^i)Z*(e'^=)e*'=(^^-^^) 

dZ{e'^^)dZ*{e'''')cosik[9i ~ 6*2]) 



(6) 



se=i/k 



The sum in the last equation can be evaluated in terms 
of the multifractal spectrum, T{q). This follows from the 
partition function approach of Halsey, et al. where 
the spectrum is defined by the implicit equation: 



^(P„)«(|AZ™|/ro)-^(^) =1 

m — 1 



(7) 



Here the surface of the cluster should be thought of as 
being divided into k boxes. The charge in box m is Pm, 
and the box size is |AZ| for that box. In our case the 
charges are all equal to 1/fc, and we must put t = — 2. 
Then we have: 



\Ak 



(8) 



where t((7) — —2 ||l^. We know that the function t{ci) 
is increasing, and that r(0) = —D. Thus 9 < and we 
confirm that the sum in Eq. (|^) is dominated by its first 



few terms, and that each scales as 



Furthermore we 



have ^\—\ ^Ylik l^fcP ~ '"o asymptotically, and the sum 
converges rapidly. 

Thus the asymptotic scaling behavior of the sum for 
the penetration length, Eq. (|4|), is the same as that of its 
terms and the sum is dominated by the first few. The 
crossover in the penetration depth must be associated 
with the crossover of the first few Laurent coefficients 
which, as we have seen, appear to behave differently from 
the prediction of Eq. (p|). The first few coefficients (see 
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Ref. Q) represent the quadrupole, octupole, etc. mo- 
ment of the cluster. We can guess that these have intrin- 
sically slow dynamics (as in the linear regime) and thus 
cross over more slowly than higher moments. Any quan- 
tity associated with averages over the charge (the growth 
zone) should share this crossover. In fact, we conjecture 
that there is one kind of crossover of lengths with the 
same correction to scaling exponent. 

We can verify this picture by numerical analysis of the 
penetration depth and other moments of the charge. We 
computed the Laurent coefficients Tq and Ai, . . . , Aiq for 
conventional DLA clusters using the method described 
earlier. For each cluster the number of points M sam- 
pling the harmonic measure was equal to the size of the 
cluster, but at least 10^; this way the numerical error in 
the computed coefficients is around 10%, and the error in 
their ensemble average is much smaller. The size of the 
clusters range from 10^ to 10^, and for each given size an 
ensemble of 1000 was taken. 

In Figure ^ we show the penetration depth and com- 
pare its crossover to that of the squared absolute value 
of the Laurent coefhcients. All of them seem to have 
asymptotic scaling of the form: 

r2(a + /3n"°-3). (9) 

The numerical results obtained by the HL method are 
close, especially for larger cluster sizes, to those shown 
on Fig. I 

What we have shown in this paper is that we can unify 
many of the puzzling results on slow crossovers in DLA. 
The conformal map approach gave us a strong indication 
that the crossovers associated with the low order Laurent 
coefficients are shared by many quantities, but that the 
asymptotic behavior is that all the characteristic length 
scales of the growth will scale as N^^^ . 

In fact, this behavior is even more general than it ap- 
pears. There is another characteristic quantity with units 
of squared length, namely the ensemble fluctuation of the 
squared radius : 

5{rl)^[<rt>-<rl>'Y/^ (10) 

where <> denotes an ensemble average. The significance 
of 5{r1) is that it is the spread, for different clusters in the 
ensemble, of r^N~'^/^ . This is not the same sort of object 
as those that we have been discussing, which are averaged 
properties of individual clusters. As was discussed in Ref. 
[Q, 5{r1) appears to scale more slowly than r^. Note, 
however, that Ref. Q used the HL method, so that only 
small N were available. 

If we go to large N we find that 5{r'^) acts in the 
same way as with the same correction to scaling: 
The apparent ensemble sharpening of the radius is also a 
crossover. This is shown in Figure |^. Once more, directly 
generating the charge distribution with random walkers 



allowed us to go to large N and reveal the crossover, 
which was not evident in the HL method. 

These results, and the last one in particular, give rise 
to the suspicion that all of the slowly scaling quantities 
in DLA growth are slaved to some underlying variable. 
Such a view was proposed some time ago by Barker and 
Ball [|l3|. In a future publication we will elaborate on 
this idea in view of the present understanding of DLA 
[Q. We think that the considerations in this paper can 
be extended to give rise to a very detailed understanding 
of DLA clusters. 
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FIG. 1. The radial displacement ry and the transverse dis- 
placement rx- The origin is at the center of charge. 
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FIG. 2. Crossover of the Laurent coefficients (open circles) 
and 5y (filled circles). 5| is approximated by the first 10 
terms of the rapidly converging sum in Eq. (Ul) . The size of 
the clusters range from 10^ to 10® in steps of VW, and for 
each given size an ensemble of 1000 was taken. 
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FIG. 3. The relative distribution width of the cluster ra- 
dius: 5{rl)/r^. The distribution width 5{rl) obeys the same 
correction to scaling as the other quantities of squared length 
dimension, see Fig. |^ The circles represent conventionally 
grown DLAs, average over 1000 clusters up to size lO**; the 
continuous lines are produced by the HL method: 400 clusters 
up to size 10'' (thick line) and 30 clusters up to size 10^ (thin 
line) . 
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